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SUMMARY 

The principal effect of subsonic trailing edges on the damping 
in roll and pitch of thin sweptback wings in a supersonic stream is 
evaluated with the aid of some conical and quasi -conical flows pre- 
viously derived. This effect is expressed in the form of approximate 
correction terms to be added to the corresponding expressions that 
are obtained when the trailing -edge disturbance is ignored. The 
results are limited to those plan forms and Mach numbers for which 
the trailing-edge disturbance does not extend beyond the leading edge 
and for which the area of mutual interference between tips and trail- 
ing edge is not large. Practical applicability Is subject to the 
limitations of linearized potential theory. 


INTRODUCTION 

At supersonic speeds the flow over a sweptback wing is closely 
related to the flow over a triangular or delta wing (fig. 1(a)) with 
the same leading-edge sweep angle. Thus, in figure 1(b), the flow 
characteristic of the delta persists out to the tip Mach cones. 

Within the tip Mach cones the flow Is altered. (See reference 1.) 

At lower supersonic speeds (fig. 1(c)), the conical Mach wave from 
the apex of the trailing edge encloses a portion of the trailing edge. 
Here also the basic delta flow is altered. (At these speeds the com- 
ponent stream velocity normal to the trailing edges Is subsonic and 
the edges are termed "subsonic.") 

The altered flow in the tip and trailing-edge regions may be 
considered to result from superposing on the basic delta-wing flow 
certain necessary additional flows. These additional flows have the 
function of canceling the delta-wing lift outside the wing tips and 
behind the trailing edge. These cancellation flows yield the dis- 
turbance pressures in the tip regions and, if the edge is subsonic, 
in the trailing-edge region. 
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According to figure 1(c), the tip and trailing-edge disturb- 
ance regions overlap . In the cross-hatched upper portion of the 
overlap, for each tip, the trailing-edge disturbance and the tip 
disturbance merely superpose. Secondary disturbances in the black 
lower portion of the overlap are caused by mutual interference 
between the tip and the trailing edge. Under fairly general cir- 
cumstances, these secondary disturbances have been found (reference l) 
to contribute a relatively unimportant amount to the integrated lift 
and moment. Accordingly, the secondary disturbances are neglected 
herein and the tip and trailing-edge disturbances are assumed to be 
independent and to superpose within the entire region of overlap. 

The problem of the tip disturbance has received a rigorous treat- 
ment in references 1 and 2 and a simplified approximate treatment in 
reference 3. The discussion herein 1 b limited to the trailing-edge 
disturbance. 

The complete flow necessary to cancel the lift behind a sub- 
sonic trailing edge is difficult to evaluate . It 1 b shown in refer- 
ence 1, however, for the case of angle of attack, that the major 
part of the disturbance pressure in the subsonic -trailing-edge region 
is obtained with only partial cancellation. (The complete cancel- 
lation is, however, carried out to close approximation in refer- 
ence 1.) This partial cancellation is afforded by a certain conical 
flow (designated flow I) in the manner shown in figure 2(a). The 
flat portion of the load distribution in section A-A represents a 
u nif orm negative lift in the triangular region behind the trailing 
edge; its magnitude is chosen to equal the value of the baBic delta- 
wing flow at the midspan. The outer portions constitute the sought - 
for pressure disturbance on the wing in the trailing-edge region. 

Similar concepts can be employed to evaluate in simple fashion 
the principal portion of the subsonic -trailing-edge pressure dis- 
turbance in the case of rolling or pitching. For rolling motion the 
basic lift to be canceled is antisymmetric. The appropriate partial 
cancellation is afforded by a certain flow HI in the manner shown 
in figure 2(b). For pitching motion the basic lift to be canceled 
increases linearly downstream. The appropriate partial cancellation 
is afforded by superposition of a constant negative lift (flow I) 
and another lift (flow IV) that increases in proportion to the dis- 
tance downstream of the apex of the trailing edge (fig. 2(c)). 

The detailed derivation of these several cancellation flows or 
"wake corrections" and others is carried through in reference 4. The 
present report is concerned with the application of some of these 
results to the determination of damping in roll and pitch for thin 
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swept wings with subsonic trailing edges. These applications take 
the form of correction terms to be added to the corresponding expres- 
sions that are obtained when the tralling-edge disturbance is ignored. 

The effect of subsonic trailing edges on the damping in roll of 
sweptback wings, which is part of the present subject, has also been 
treated recently in reference 2. The treatment therein achieves 
nearly exact cancellation of the lift behind the trailing edge by the 
superposition of infinitely many infinitesimal flows. The procedure 
is a development of the approach used in reference 1. The present 
treatment differs in employing approximate cancellation by super- 
position of a single flow. This simplification permits giving the 
over-all result compactly in closed form. 

The results presented herein are applicable for those super- 
sonic speeds for which the trailing-edge disturbance does not extend 
beyond the leading edge. At the lower Mach numbers for which the 
leading edge is enveloped further correction is required. This low 
Mach number problem is studied for the case of an gle of attack in 
reference 5. 


CORRECTION FOR DAMPING IN ROLL 


For a sweptback wing with subsonic leading edges rolling with 
angular velocity p, the basic delta-wing lift distribution to 
be canceled behind the trail i ng edge is (equation (26), reference 6) 


2 

2(x+c_.)ypm I(m) 

C P,1 - } ' , 

0V (x+c r ) m 2 -P 2 y 2 


where 

l(m) = 2(1^) 

(2-m 2 )E«(m)-m 2 K» (m) 

and x is measured from the apex of the trailing edge. (All symbols 
are defined in appendix A.) This load distribution has the anti- 
symmetric character shown in figure 2(b). For the inexact cancel- 
lation employed herein, the distribution is approximated by its 
tangent at y=0: 
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2>C 


C . ~ 


P>A 


P>A y dy 


'y=0 


^ 2pml(m) 

0+ 

PV 


( 2 ) 


(The limitations on the use of this approximation are considered in 
the DISCUSSION . ) 

The load-distribution equation (2) can be canceled by flow III 
of reference 4, which can be written as 


C p = Cy 


(3) 


where 


C = constant 


in the region 0 <. | c| S. n behind the trailing edge of the wing. 
Thus the choice 


2pml(m) 

PV 


(4) 


achieves the desired cancellation. 


Cancellation flow TTT includes also a lift disturbance ahead of 
the trailing edge in the region n <. |Ol <.1. Beference 4 gives for 
this disturbance, in the present notation. 


C = Cx 

p 


a [ic^kJ-nS-C^hj] - A^(l-Cf 2 )(a 2 -n 2 ) 
P [i(i:)-n 2 K(l0l 


( 5 ) 


where 

C = same constant as in equation (3) 


4 = sin" 1 A 1^- 


1-n 
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k = *\l 1-n 2 

= Py/x 

The subsonic -trailing -edge correction to the damping in roll results 
from integrations to obtain, the contribution of equation (5) to the 
rolling moment of the wing. 

The correction is evaluated as follows: Equation (5) is written 

in the form 


Cp = Cxf(c) 

Then the incremental rolling moment on an elementary triangular region 
(hatched area’, fig. 3) is given by 


d 



pb 
n 2 a 

x=0 


yCxf(a) i di ^ 

p 


cpVrCa) da 

ns - 

64 O 5 


The incremental rolling moment experienced by the finite 
triangular region n < a <L 1 of figure 3 is given by 



2 4 
CP b 

64 


^1 





( 6 ) 


The function f (a) is the expression in braces in equation (5). 
Thus the Integrand in equation (6) involves elliptic integrals and 
related functions. The integration may be effected without diffi- 
culty with the aid of the elliptic -function substitutions of 
appendix B. The result is 

AL' = _ Cpb 4 
Tj- pV 2 64n 


|_ 4 ECkJ-n^OOJ 


( 7 ) 
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The following standard definitions are convenient: 

E(k) = E'(n) 

E(k) = Z* (n) 

With these substitutions, equation (7) becomes 

Al* Cpb 4 

| pY d 64n 


II - - 


1-n 


4 E'(n)-n 2 K« (n) 


( 8 ) 


The value of C given by equation (4) is to be substituted back 
in equation (8). Then the expression may be converted to coefficient 
form by division by Sb and by pb/27. Multiplication by a factor 
of 2 accounts for both wing panels on the complete wing. The final 
result is (with both sides multiplied by 3) 


4«B, ) - 

p 8n 


* 1-n 2 

4 E’ (n)-n 2 K l (n) 


(9) 


in which the loading- and ‘trailing-edge sweep parameters m and n, 
respectively, are related by 


1 a 1 _ 1-A 

n ~ m ' ap I+X 


do) 


where A is the aspect ratio and A the taper ratio of the wing. 

Equation (9) gives a correction far the damping-in -roll deriv- 
ative C 7 when the trailing edge of the wing is subsonic. This 
P 

correction is algebraically additive to the value of PC, computed, 

V 

as in reference 3 (in effect), by integration of the basic delta- 
wing lift distribution together with the primary tip correction. 


Applicability of equation (9) is limited to those plan forms and 
supersonic speeds for which the trai ling-edge disturbance does not 
envelope the leading edge of the wing. This limitation is a conse- 
quence of the shape of the region of integration assumed in fig- 
ure 3. The limitation may be expressed by 
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or alternatively "by 
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— ~ < - < 1 + 
A ft 1 4- A m " 


» A / 1 . % \ 


(U-t) 


CORRECTION FOR DAMPING IN PITCH 

In the case of a sveptback wing with subsonic leading edge 
pitching about its apex, the basic delta-wing lift to be canceled 
behind the trail i ng edge is (equation (29), reference 6): 



4qG(m) [jJfx+Cp ) 2 m 2 -p 2 y 2 [] 

pv 


( 12 ) 


where 


G(m) 


1 C 

-m 


(l-2m 2 )E’ (mJ-Hn^K 1 (m) 


and x is measured from the apex of the trailing edge. Now the 
variation of C p ^ with y is very gradual for a certain range 

of yj for the approximate cancellation employed herein, this vari- 
ation is neglected. Thus the lift canceled is given by equation (12) 
with y therein set equal to zero. 



8c r qmG(m) 8qmG(m) 

o + x 

PV PV 


(13) 


(The replacement of the exact curve (equation (12)) by its tangent 
(equation (13)) should give unusually accurate results because the 
first three derivatives of the function with respect to y vanish 
at y =* 0. The limitations on the use of this approximation are con- 
sidered in the DISCUSSION.) By suitable choice of constants, the 
first term can be canceled by flow I of reference 3 'and the second 
term by flow IV of reference 4, which are of the form 
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flow Is C p = C x 
flow IT: C = Cj-yi 


0 < a < n 


Thais the choice 


"IV 


8c r qmG(m) 

PV 


8qmG(m) 

PV 


(14) 


(15) 


achieves the desired cancellation* 


Equations (14) refer to the portions of the cancellation flows 
hehind the wing trailing edge (that is, 0 < |a|< n). Each of these 
flows also includes a lift disturbance ahead of the trailing edge in 
the region n < |<j! < 1. Beference 4 giveB for these disturbances 


flow Is C a C T 

P I K(k) 


flow IV 


. c = c „x 

p vr [_ K(k)-E(t) J 


> n < o < 1 (16) 


The subsonic-trailing-edge correction to the damping in pitch results 
from integrations to obtain the contribution of equations (16) to the 
pitching moment of the wing. 


The pitching-moment correction referred to the apex of the wing 
is desired, but it is more convenient at first to refer the moments 
to the apex of the trailing edge (y-axis of fig. 3). This procedure 
entails obtaining both lifts and moments, inasmuch as the lift is 
involved in the transfer of the moment reference axis back to the 
wing apex. The first of equations (16)' (flow I) is of the 
form C p = Cjf(a). The lift and the moment imparted by flow I to 

the elementary hatched triangle of figure 3 are therefore 
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Pb 

r 2a 

Cjf (a) x dx ~ 

4 


Cjpb f(a) da 
8 a 2 


d 



Pb 

~ 2a 




xCjf (a) x dx 


da 

P 


Cj-ph^io) da 
24a 3 


The lift and the moment experienced by the finite triangular area of 
figure 3 are obtained by integrating between the limits n and 1. 
The function f (a) therein is the coefficient of C-j. in equa- 
tions (16) o The substitutions of appendix B are helpful in effect- 
ing the integrations* The results are 



(17) 


(18) 


The second of equations (16) (flow IT) is of the form 
Cp = Cjyxf (a)* The lift and the moment imparted by flow IT to the 

elementary hatched triangle of figure 3 are therefore 
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Pb 

~2a 


Cryzf(a) i cLx 


C-j^pVf (a) da 


24a 


d 



Pb 

r 2 a 




XCjyXf (a) I dx -p- 


-c^pV^a) da 


Tha lift and the moment experienced by the finite triangular area of 
figure 3 are obtained, as before, by integrating between the li m its n 
and 1. The function f(a) is the expression in brackets in the 
second of equations (16). The result of the integrations is 



c^pV 

48n 2 


E'(n)-n 2 g' (n) 
K* (n)-E’ (n) 


A IV M 




n 1-n 2 
4 K* (n)-E’ (n) 


( 19 ) 


( 20 ) 


All the data are now at hand for calculating the incremental 
pitching moment referred to the wing apex. First the incremental lift 
will be needed. 


AL = A I L + A IT L 
| PY 2 | PY 2 | PY 2 


( 21 ) 
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Then the incremental pitching moment is 


H 
Oi 

o 
o 

The several moments and lifts in these two equations my "be 
taken from equations (17) to (20 ) , using for and C^. therein 

the values determined in equations (15). Then the left sides may he 
converted to coefficient form hy dividing equation (21) by S 
and qc/2V end equation (22) by Sc and qc/27. Multiplication by 
a factor of 2 accounts for both wing panels on the complete wing. 
The final results are, with both sides multiplied by 6, 


AM 

1 o 

2 PV 2 


A-j-M 


AjyM 


' °r 


AL 


( 22 ) 


2 P Y 2 


2 


tv) 


\ ApmG(m) j ^ c r 

1 n/2 

2 CAP 

1 E'(n)-n%'(n) 


i)” » f • 

(n) 

3 c n 

K' (n)-E'(n) 



(23) 


k K) 


APmG(m) 

r 4 ^ 

1 11/2 

n 

L 

K* (n) _ 


2 c Ap 

3 c c n 

2 2 2 
1 c A P 

6 ? n 2 


2 - 


E«(n) 

K* (n) 


E'tnJ-n^' (n) | 
K’ (n)-E' (n) | 


11 - - 


1-n 


4 E« (n)-E* (n) 


(24) 


Equation (24) gives the value of a correction for the damping- 
in-pitch derivative C_ when the trailing edge of the wing is sub- 

sonic. This correction is algebraically additive to the value 
of C„ computed by integration of the basic delta -wing lift distri- 


m_ 


button alone. Equation (23) gives a corresponding corr.ecticn for the 
less -important lift-due-to-pitc hin g derivative Cj^. Both equations 

are subject to the Mach number-plan form limitations expressed by 
equations (11a) and (lib). 
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Equations (23) and (24) refer to a wing pitching about its 
apex with moments referred to the apex. The following transformations 
will provide values of the corrections to the derivatives (designated 
by primes) referred to an axis a distance h downstream of the apex: 



The quantities A(BC r \ and A^BC m ^ on the right side are given by 


t\) TV , , 

equations (23) and (24). The quantities A^BC^ j 
be obtained from reference 1 (corrected in the errata sheet thereto). 


and A 


K) 


may 


Similar procedures for the treatment of the effect of subsonic 
trailing edges on the derivatives Cj (effective dihedral) and Cj 


are considered briefly in appendix C. 


DISCUSSION 

The circumstances under which the approximations used herein lead 
to acceptably small error in the calculation of the wing load in g 
require examination. There are two points of departure from " exact” 
linearized potential -flow theory employing the Kutta condition at the 
trailing edge: (1) the basic delta -wing load distribution behind the 

wing is canceled only approximately by superposition in each case of 
a special flow, and there remains a generally small residual lift in 
the wake (fig. 2 ); (2) the mutual interaction between the tips and 
the trailing edge (black regions of fig. l) is neglected. 

In more elaborate treatments (reference 1 for angle of attack, 
reference 2 for rolling motion) the lift in the wake is canceled 
almost completely, and the mutual interaction between the tip and the 
trailing edge is also taken into account. In reference 1 two examples 
are presented from which the separate contributions of the several 
degrees of approximation may be assessed. (Similar examples in refer- 
ence 2 are not so well Buited to such an assessment, because of an 
imperfect correspondence. One of these examples is, however, here- 
inafter considered.) The examples refer to an un tapered wing of 
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aspect ratio 1.72 and a tapered wing of aspect ratio 3.85 and taper 
ratio 0.179, 'both with 63° leading-edge sweepback. The Mach number 
is 1.5. The breakdown for lift is as follows: 


Source of Lift 

TEa tapered 
wing 

(percent) 

Tapered 

wing 

(percent) 

(l) Basic delta -wing loading 

100.0 

100.0 

(2) Tip effect 

-14.7 

-1.6 

(3) Symmetrical wake correction 

-12.1 

-2.7 

(4) Oblique elements in wake 

-.82 

-.59 

(5) Tip -trailing -edge mutual interference 

3.0 

.04 

Total lift 

75.4 

95.1 


The "symmetrical wake correction" is cancellation flow I of this 
report (fig. 2). The "oblique elements in the wake" is an addi- 
tional flow that, added to flow I, (almost exactly) completes the 
cancellation of the basic delta-wing lift behind the wing. Thus for 
the untapered wing the item of -0.82 percent represents the error 
in computed wing loading associated with the failure of flow I to 
cancel completely the basic delta -wing loading behind the wing. 
Finally, the item 3.0 percent represents the error associated with 
neglect of the mutual interference between the tips and the trailing 
edge. The errors of the present treatment therefore amount to 
-0.82 + 3.0 or 2.2 percent for the calculation of the lift of the 
untapered wing. (Note that these errors are expressed as percentages 
of item 1.) 


The examples afford a similar breakdown for the pitching moment 
about the wing apex: 


Source of Moment 

Tkrtapered 

wing 

(percent) 

Tapered 

wing 

(percent) 

(l) Basic delta-wing loading 

100.0 

1 

(2) Tip effect 

-23.9 


(3) Symmetrical wake correction 

-17.7 


(4) Oblique elements in wake 

-1.37 


(5) Tip-tralll ng-edge mutual interference 

5.3 

■ 

Total pitching moment 

62.3 

. 92.3 
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The errors of the present treatment for the calculation of the pitch- 
ing moment for the untapered wing therefore come to -1.37 + 5.3 or 
3.9 percent and for the tapered wing to -0.95 percent (again "based on 
item (1) = 100 percent). 

It is possible to infer, in a general way, from these examples 
the circumstances under which the approximations of this report, 
namely, the neglect of items (4) and (5), are reasonably satisfactory. 
The difference between cancellation flow I and the basic delta-wing 
lift in figure 2 corresponds to item (4). For a given section A-A 
across the wake, this difference decreases with the ratio QR/PS. 

For the evaluation of item (4) section A-A may be taken to pass through 
the trailing ends of the wing tips (that is, QR equals wing span b 
and P and S lie on the leading edges extended). The ratio QR/PS 
for this section is designated T. In the untapered-wing 
example, T = 0.629; in the tapered -wing example, T = 0.926. 

Now consider variations of wing sweep, taper, and Mach number 
that leave T constant. It is expected that item (4) will change 
approximately in proportion to item (3) (exactly, if the trailing- 
edge-sweep - Mach number parameter n is unchanged) . On this basis 
.item (4) in the examples may be adjusted to equal 1 percent and 
item (3) scaled up accordingly. The results are 



Item (3) in 

T 

C L 

Crr, 

m 


a 

a 


(percent) 

(percent) 

0.629 

-14.7 

-12.9 

.926 

-4.6 

-4.2 


These four points are plotted in figure 4. A straight -line variation 
of T against item (3) is not inconsistent with the points. The line 
drawn in is considered to represent approximately the condition that 
item (4) equals 1 percent. Accordingly, for values of T below the 
line the neglect of item (4) will yield an error of less than 1 per- 
cent (based on basic delta -wing loading = 100 percent). The par- 
ameter T may be calculated from the formula 


1 



— + — — — — 
m AB(l+?0 


( 26 ) 


1300 



NASA TN 2146 


15 


A corresponding criterion for the neglect of item (5), the mutual 
interference between tips and trailing edge, c ann ot he formulated with 
much confidence. The ratio of the area moment of the wing regions 
affected by this interference (black regions in fig. 1) to the area 
moment of the entire wing surface is presumed to be a rough guide in 
the estimation of moments. In the tapered -wing example both item (5) 
and this ratio for the area pitching moments are negligible. In the 
untapered -wing example this ratio for the area pitching moments 1b 
17 percent and item (5) is 5.3 percent. For the same wing this ratio 
for the area roll in g moments is 15.4 percent and item (5) for ro llin g 
motion is 4.3 percent. (See subsequent example, taken from refer- 
ence 2.) If the correspondence can be stretched somewhat, the rule 
that the percent error in calculated moment associated with neglect of 
item (5) is about one -third the percent area -moment ratio may be used 
as a tentative criterion. Equations (11a) and (lib) constitute an 
additional limitation that states, in effect, that the black inter - 
ference regions in figure 1 may not impinge on portions of the leading 
edge. This limitation was introduced partly for convenience in the 
integrations; both this and the area -moment limitation act in the Bame 
direction. Of the two, that one governs that results in the smaller 
area -moment ratio. 


These several criterions for the suitability of the approximations 
of this report were inferred from examples worked out for a wing at an 
angle of attack. The similarity of the features of the load cancel- 
lation for rolling motion and for pitching motion (fig. 2) suggests 
that these criterions are of the right order for these more general 
motions. (Note that the appropriate axis, in each case, should be used 
for the area -moment criterion.) An indication that the criterion 
for T is not unconservative for the case of rolling comes from the 
following calculation of for the untapered wing of reference 1, 

taken (except for item (3)) from reference 2: 


Source of Moment 


Un tapered 
wing 


(l) Basic delta-wing loading (reference 2) 

‘ (2) Tip effect (reference 2) 

(3) Partial trailing-edge correction (equation (9) herein) 
(3a) Total trailing-edge correction (reference 2) 

(5) Tip-trailing edge mutual interference 


(percent) 

100.00 

-37.54 

-8.06 

-8.23 

4.26 


Item (3) is the approximate trailing-edge correction computed by 
equation (9) of the present paper. Item (3a) is the trailing-edge 
correction computed in reference 2 by substantially exact cancellation 
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of the residual .lift in the wake . Thus the difference "between 
items (3a) a n d (3) corresponds to item (4) of the earlier examples. 
This difference is very small: only 0.17 percent of the uncorrected 

value of Cj , item (1). This value for item (4) is considered to 
P • 

have some uncertainty in view of the fact that item (3a) was obtained 
"by numerical integration. If the point T = 0.629, 
item (3) = 8.06 percent, is plotted on figure 4 it falls, as it 
should, within the region for which item (4) is less than 1 percent. 

So far the approximations in this report have "been Justified 
solely on the "basis of their closeness to the values required "by the 
"exact" linearized potential theory. For real wings, however, the 
action of viscosity in producing a "boundary layer partly vitiates 
the predictions of potential theory. The houndary-layer flow for 
sweptback wings is such that the deviation from potential theory is 
expected to he appreciable in the vicinity of the wing tips and the 
trail ing edge. Thus, the approximations that have been introduced' 
herein are probably within the accuracy of applicability of the 
potential theory itself, so long as the several criterions have been 
met. 


Lewis Flight Propulsion Laboratory, 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, January 6, 1950. 
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APPSSDH A 


SYMBOLS 

The following symbols are used in tills report: 
A aspect ratio 

h span 

C,Cj,Ciy constants 

C T coefficient of lift due to pitching 


a 


m 


lift 


1 py 2 s to 

2 H VST / 


lift-curve slope 


lift 
| PV 2 Sa 


coefficient of damping in roll 


rolling moment 


1 2 / Pt> 


pb\ 

27/ 


coefficient of rolling moment due to yaving 
rolling moment 


I p v 2 si> (§) 


coefficient of rolling moment due to sideslip 
/ rolling moment N 


l 


| p7 2 ShP 


coefficient of damping in pitch 


pitching moment 


1 pt 2 sc ( SZ) 

2 V, 27/ 
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“a 


°r 

35(0,k) 

E(k) 

E’ (m) 

E«(n) 

E(0,k) 

f(o) 

G(m) 

h 


coefficient of pitching moment due to angle of attack 
^ pitching moment^ 

1 2 
I PY Sea 

pressure coefficient (lift per unit area divided "by 
| PV^Sca) ' 

rW 2 


mean aerodynamic chord 

““ rV 2 


l 

s 


(local chord ) c dy 


average chord 
root chord 


1 

h 


-h/2 


-h/2 

(local chord) dy 


incomplete el li ptic integral of second kind with 
amplitude 0 and modulus k 

complete elliptic integral of second kind with 

'modulus k 


E 


E 


(v^) 

(v^) 


- E (f> k ) 


incomplete elliptic integral of first kind with 
amplitude 0 and modulus k 

function of <7 defined in text 


1-m 


(l-2m 2 )E' (m)+m 2 E» (m) 

distance of center of gravity downstream of wing apex 


a 
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I(m) 

2(l-m 2 ) 

(2-m 2 )E' (mJ-m^K* (m) 

K(k) 

« 

complete elliptic integral < 


modulus k = F^,k^j 

K* (m) = 

K (V 1-m 2 ) 

K* (n) = 

K ( V 1-n 2 ) 


k 

L 


modulus of elliptic integral 
lift 


(= f\J 1-n 2 ) 


L« 

M 

m 

n 

P 

Q. 

r 

8 


rolling moment (positive in sense of right-hand screv 
proceeding upstream) 

pitching moment (positive in sense of right-hand screw 
proceeding spahwise to right) 

P times contangent of leading-edge sweep angle measured 
from y-axis 

P times contangent of trailing-edge sweep angle measured 
from y-axis 

angular velocity of rolling 
angular velocity of pitching 
angular velocity of yawing 
wing area 


u,v,w disturbance velocity components along x-, y~, and z-axes, 
respectively 

V stream velocity 

Cartesian coordinates: x-axis parallel to free -stream 

direction; y-axis horizontal and toward the right, 
looking upstream; z-axis vertically upward 


x,y,z 
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APPENDIX B 


AIDS TO INTEGRATION 


The integrations in this paper that involve elliptic integrals 
are most easily evaluated with the aid of the following elliptic 
f unction substitutions, together with the tables of integrals in 
references 7 and 8, and integration by parts: 

F(0,k) « u 

pu 

E(0,3c) = dn^u du 

vjo 

^1-a 2 «= k sn u 
t\jo 2 -n 2 k cn u 


where 


a *» dn u 

dCJ = -k 2 sn u cn u du 


0 « sin"- 1 - 


k = Nl~n 2 
k* » Vl-k 2 = n 
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APPENDIX C 


EXTENSION TO C, AND C 7 
P r 


The antisymmetric component of the basic delta-wing loading 
is qual i tatively similar for rolling p, sideslip 6, and 
yawing r. Thus the procedure developed herein for the damp ing- 


in -roll derivative Cj may he adapted without difficulty to the 

effective-dihedral derivative C 7 and. once the correct delta - 

ip 


wing loading is known for yawing, to the derivative C 7 • This pro- 

*r 

cedure will, however, he less accurate for C 7 than for C 7 , 


because the delta-wing loading to he canceled in the wake is less 
well approximated by flow IH. An approximation that is closer for 
sideslip, hut still simple, is a flow specified by u = Cy/Ccy+x) 
in the region 0<la|< n and w = 0 in the regions n<la| <1. 
This flow is difficult to solve far u in the w = 0 regions, hut 
a flow closely similar in the right- hand w = 0 region may he 
readily solved by a procedure in reference 4. This similar flow is 
obtained by replacing w = 0 by u = 0 in the left-hand 
region -l<o<'-n. The f in al solution for u in the right-hand 
region n< CT <1 is to be applied for both nS.<J <1 and -1 < 0 <-n - 



1300 


NACA TN 2146 


23 


REFERENCES 

1. Cohen Doris: The Theoretical Lift of Flat Svept -Back Wings at 

Supersonic Speeds. MCA TN 1555, 1948. 

2. W a l ke r, Harold J., and Ballantyne, Mary B. : Pressure Distribution 

and Damping in Steady Roll at Supersonic Mach Numbers of Flat 
Swept-Back Wings with Subsonic Edges. MCA TN 2047, 1950. 

3. Waive s tut o, Frank S., Jr., Mar golds, Kenneth, and Ribner, 

Herbert S. : Theoretical Lift and Damping in Roll of Thin 

Sweptback Wings of Arbitrary Taper and Sweep at Supersonic 
Speeds. Subsonic Leading Edges and Supersonic Tra iling 
EdgeB. MCA TN 1860, 1946. 

4. Ribner, Herbert S. : Some Conical and Quasi -Conical Flows in 

Linearized Supersonic Wing Theory. MCA TN 2147, 1950. 

5. Cohen, Doris: Theoretical Loading at Supersonic Speeds of 

Flat Swept-Back Wings with Interacting Trailing and 
Leading Edges. NACA TN 1991, 1949. 

6. Brown, Clinton E. , and Adams, Mac C. : Damping in Pitch «nd Roll 

of Triangular Wings at Supersonic Speeds.. MCA Rep. 892, 1948. 

7. Peirce, B. 0.: A Short Table of Integrals. Ginn and Co., 3d ed., 

1929, p. 18. 0 

8. Jahnke, Eugene, and Emde, Fritz: Tables of Functions with Formulas 

and Curves. Reprint of 1943 Edition, Dover Pub. (New York), 
pp. 96-97. 



CO 




(a) Delta wing. (b) Sweptback wing formed 

from delta wing. 



(c) Sweptback wing formed 
from delta wing at lower 
supersonic Mach number. 


Figure 1. - Related wings in supersonic flow. Hatched regions explained in text. 
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Figure 2. - Sweptbaclc wing showing basic delta -wing load 

distribution along section A-A and approximate cancellation 
of this load behind wing by superposition of special flows. 
Special flows are plotted with reversed sign. Cross- 
hatched areas represent Induced changes in loading ahead 
of trailing edge. 
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Figure 4. - Estimated variation of t with trailing-edge correction 
computed by present method under condition that error in correc- 
tion, exclusive of tip interaction, equals 1 percent of contribu- 
tion of basic delta-wing loading. 
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